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Abstract 

A Large Ion Collider Experiment (ALICE) is one of the four large experiments at the Large Hadron Collider 
(LHC), and it is dedicated to the study of ultra-relativistic heavy-ion collisions, with the goal of investigating the 
properties of the high-density state of QCD matter produced in these collisions. 

The study of D meson production azimuthal anisotropy and the measurement of their elliptic flow (v'2) can provide 
insight on the degree of thermalisation of charm quarks in the medium and on the charm hadronization mechanism. 

We present the measurement of the D + and D° meson V2 in Pb-Pb collisions at ysNN = 2.76 TeV at the LHC with 
ALICE. We discuss the details of the analysis and we show the results obtained from data samples collected in 2011. 
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1. Introduction 

At high energy densities, Quantum-Chro mo-Dynamics predicts the formation of a state of matter in which partons 
are not confined anymore inside hadrons. This state of matter is called Quark Gluon Plasma (QGP). A Large Ion 
Collider Experiment (ALICE) HI is one of the four large experiments at the Large Hadron Collider (LHC) and its 
main experimental goal is the study and characterisation of the QGP created in the high energy Pb-Pb collisions at 
ysNN = 2.76 TeV delivered by the LHC. In ALICE, vertexing and charged particle tracking at central rapidity relies 
on six layers of high resolution silicon detectors of the Inner Tracking System (ITS), with the two innermost layers 
equipped with Silicon Pixel Detecetors (SPD), a large volume Time Projection Chamber (TPC), and a high granularity 
Transition-Radiation Detector (TRD). In non-central Pb-Pb collisions the fireball shows azimuthal anisotropy with 
respect to the reaction plane (the plane defined by the beam direction and the impact parameter). Due to collective 
effects, this spatial anisotropy is translated into a momentum anisotropy driven by different pressure gradients in the 
in-plane and out-of-plane regions. This anisotropy can be quantified expanding the particles momentum azimuthal 
distribution (with respect to the event plane) in a Fourier series. At mid-rapidity the main coefficient is the second 
coefficient of the Fourier expansion v'2, also called elliptic flow [2|. The measurement of heavy-flavour particles V2 
is particularly interesting as it can probe the degree of thermalisation of the medium, the interaction of heavy quarks 
with the QGP and the heavy flavours hadronization mechanism J3]lll. The charm v'2 has been measured at RHIC by 
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the PHENIX collaboration via non-photonic electrons [5|. Thanks to its vertex reconstruction, tracking and particle 
identification capabilities, ALICE can perform a direct charm V2 measurement through fully reconstructed D meson 
hadronic decays. In these proceedings we present the result on the measurement of D meson V2 from ~ 9.5 ■ 10 6 
(~ 7.1 ■ 10 6 ) Pb-Pb collisions in the centrality class 30-50% (15-30%). 



2. Event selection and event plane determination 

The data sample used for the D meson analysis has been collected with a minimum bias trigger given by the 
coincidence of signal in the VZERO scintillators and the SPD, and a semi-central trigger tuned to select events with 
centrality < 50% (where centrality is associated to a central event). The centrality selection of the events is done on 
the distribution of signal in the VZERO scintillators. A Glauber model fit is performed on MonteCarlo simulations to 
extract the total cross section and normalise [6 |. The event plane is measured from the distribution of charged tracks 
in the < rj < 0.8 region of the ALICE TPC using the equation 



1 

¥ = - tan" 
2 
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(1) 



where *P is the second harmonic event plane and 0; is the angle of the 2 th track in the ALICE reference frame 0. 
Weights are applied to the tracks depending on their angle in the ALICE reference frame to account for different 
efficiencies among TPC sectors. To compute the event plane resolution each event is splitted in 2 sub-events with 
each track in the < rj < 0.8 range randomly assigned to one of the sub-events and the sub-event plane angle ^a/b 
is computed. From the two sub-events plane angles the resolution is computed following the prescriptions in |7|. 
The resulting event plane is flat to about 2% level and its resolution is 0.86 in the centrality class 30-50% and 0.9 in 
15-30% centrality. 



3. D meson reconstruction and V2 extraction 

Open charm hadrons reconstruction at ALICE is based on the invariant mass analysis of fully reconstructed 
hadronic decay topologies. The D meson vo has been measured for the decay channels D + — > K~7t + tt + , D°-> JT TV 
and the D* + — > D°7r + analysis is ongoing. As the mean proper decay length ct of these particles is of the order 
of ~ 150 (D°)— 300 (D + )/im, the strategy to reduce the large combinatorial background is based on the selection of 
displaced-vertex topologies, i.e. with large separation between the primary and the secondary vertices and on the 
good alignment between the reconstructed D meson momentum and the flight line. Particle identification is provided 
by time of flight measurements in the Time Of Flight (TOF) detector and by specific energy deposit in the TPC. The 
PID selection strategy was defined with the aim of preserving all the signal. Further analysis details can be found in 
0. The invariant mass distribution of the reconstructed D meson candidates is splitted in the in-plane and out-of- 
plane regions defined as the regions |o < A(p < f ) U [t < A<p < 7rj and [ j < A<p < respectively, where A<p is the 
azimuthal angle of the reconstructed candidate with respect to the event plane. An invariant mass analysis is used to 
extract the amount of signal candidates measured in the two regions (Njn, A^out)- A Gaussian distribution is assumed 
for the signal shape. In each pj bin the width of the Gaussian is fixed to the value obtained fitting the A(f> integrated 
invariant mass distribution. From the values of Nin and jVout it is possible to directly compute the elliptic flow as 

n N !N - Nout 

V? = . (2) 

4Ni N +Nout 

Other methods to extract v>2 based on 2-particles correlations |9| have been implemented. Their results have been 
found in agreement with the v>2 obtained from equation [2] 

4. Results 

4.1. Systematic uncertainties 

Several sources of systematic uncertainties have been considered in the analysis. The main contributions to the 
total systematic uncertainty come from uncertainties on the yield extraction and from the topological cut selection. 
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Figure 1: D° (blue circles) and D + (grey crosses) V2 
measured by ALICE in the 30-50% centrality class. 
Empty boxes show systematic uncertainties from data, 
shaded area the uncertainty on the B feed-down. Black 
points are charged hadron V2. 
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Figure 2: D V2 measured by ALICE in the 30- 
50% (blue circles) and 15-30% (red squares) centrality 
classes. Empty boxes show the systematic uncertain- 
ties from data, shaded area the uncertainty on the B 
feed-down. 



The systematic uncertainty on the yield extraction is estimated using different background functions in the fitting 
procedure and allowing the Gaussian width to vary between in-plane and out-plane invariant mass distributions. The 
analysis was also repeated using three different sets of topological cuts. Yield extraction and cut variation gave a total 
systematic uncertainty of absolute value 0.02-0.05 each (depending on the channel and on px) for the D° analysis. 
Other minor sources of systematic uncertainties are the centrality dependence of the resolution in the centrality range 
considered, that introduces a systematic uncertainty of ~3% and a 7% asymmetric systematic uncertainty towards 
larger values of V2 that comes from the definition of the sub-events. 

4.2. B feed-down subtraction 

The signal sample contains a fraction of D meson coming from B decays, thus the measured elliptic flow v* 8 
is a combination of the elliptic flows of prompt D (v^ romp ) and D from B feed-down (yf eed - down ) mesons that can be 
expressed as — / pl0 mpt v plompt + ( 1 - /prompt )v2 eed ~ down , wnere /prompt is the fraction of prompt D meson in our sample. 
From MonteCarlo simulations and FONLL 1 1 1 predictions we estimated /prompt to be in the range 0.7-0.95 depending 
on channel and p T , as well as on the topological cuts applied. The relative ^aa(a) - <t aa > ^r^~ suppression of 
feed-down and prompt candidates also affects the value of /prompt- The effect was estimated by computing the prompt 
fraction varying this ratio in the range 0.5 < rP™ 11 ? 1 / r^™" < 2 Most of the mo dels predict < v ^ ed - down < v prompt 
due to the larger mass of the b quark. For our measurement we make the assumption that yf eed - down = yP rompt . This 
assumption introduces an asymmetric systematic uncertainty on our measurement towards higher V2 values. The 
maximum value of this uncertainty is given by vP rompt // P rompt, that corresponds to the limit in which yf eed down = 0. 

4.3. D meson elliptic flow 

The D meson elliptic flow was measured at ALICE for the D° — > K ~n + and D + — > K~n + n + channels in the 
centrality class 30-50%. The results are shown in figure [T] for several pj bins in the range 2 < pj < 16 GeV/c. In 
the range where also D + measurement is available (3 < pi < 8 GeV/c), the two measurements are compatible within 
statistical uncertainties and both are compatible with the ALICE measurement of charged hadron V2 IfTTI in the same 
rapidity region. The D° V2 measurement was also carried out in the centrality class 15-30% for 2 < pi < 16 GeV/c. 
The result is shown in figure [2] compared to the result in 30-50%. Our measurement has been compared to theoretical 
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predictions fl21[l3][T4][T5]], shown in figure [3] (left). The same models have been compared with the D meson Raa 
measurement from ALICE (figure [3] right). 
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Figure 3: Left: D° v% measured by ALICE in the 30-50% centrality class (blue circles). Empty boxes show the 
systematic uncertainties from data, shaded area the uncertainty on the B feed-down contribution. Right: D meson Raa 
measurement O (black). The two measurements are compared to theoretical models lfl2l [T3l fT4l [TBI . 



5. Conclusions 

The measurement of D meson V2 at ALICE from Pb-Pb collisions at y^NN - 2.76 TeV have been presented. The 
results indicate V2 > (with 3<x significance) in the pj range 2 < pj < 6 GeV/c and are consistent among different 
mesons and different methods to compute \i- The result is also compatible, within uncertainties, with the V2 of charged 
hadrons in the same p T and rapidity region. An indication of dependence of v-i with centrality was observed between 
the centrality classes 15-30% and 30-50%. Models of heavy quark transport in the medium can describe the data. 
However, a simultaneous description of D meson \>2 and Raa is still lacking. 
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